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Abstract 

The DO collaboration has recently introduced new variables, a-r and </>* to more accurately 
probe the low Qt domain of Z boson production at hadron colliders than had been pre- 
viously possible through a direct study of the Qt distribution. The comparison of such 
accurate data to precise theoretical predictions from QCD perturbation theory will yield 
important information on the ability of resummed QCD predictions as well as parton shower 
models to describe the low Qt domain and should enable more stringent constraints on non- 
perturbative effects. In the present Letter we provide analytical predictions for the above 
mentioned variables, that contain resummation of large logarithms, including next-to~next- 
to leading logarithmic (NNLL) terms, supplemented by exact 0(a;^) next-to-leading order 
calculations from MCFM. 



1. Introduction 



The transverse momentum distribution of lepton pairs produced via the Dreh-Yan process jH at 
hadron colhders, has been one of the most extensively studied variables in QCD fl-Q. In spite of 
a wealth of literature on the subject, measurements of Qt spectra of lepton pairs (or equivalently of 
gauge bosons that decay to lepton pairs) and various theoretical approaches to the low Qt regime still 
form a significant part of modern QCD phenomenology (for recent theoretical and phenomenological 



work on the subject see for instance Refs. [l5|, llSl. Il7l|) 



The reason for strong and sustained theoretical interest in the low Qt domain of say Z production 
at hadron colliders is because this region receives important contributions from the multiple emission 
of soft and/or collinear gluons. The ability of the theory to make precise predictions for the low Qt 
domain, which incorporate the resummation of large logarithms, can thus be taken as a strong signal 
that the relevant perturbative QCD dynamics is well understood. The low Qt region is hence also of 
interest for the purposes of testing, further developing or improving parton shower models embedded in 
QCD event generators and a comparison between parton showers and analytical resummed predictions 
is always of interest in this context. 

Moreover, the state of the art for perturbative calculations in QCD is currently such that for variables 
such as the Drell-Yan and Higgs boson Qt spectra resummed predictions exist to NNLL accuracy in the 
resummed exponent, which when combined with NLO estimates from fixed-order calculations ought to 
produce some of the most precise predictions for QCD observables to date. 

However, while the existence of such accurate resummed calculations is encouraging, there are a 
number of issues that are omitted from the perturbative resummation which could play a significant role 
in the confrontation of theory with experiment. One aspect of physics at hadron colliders, that becomes 
important at small Qt, is the role of non-perturbative effects commonly attributed to the intrinsic fc^ 
of partons within the proton. One may therefore view any opportunity to compare precise perturbative 
predictions with accurate experimental data as a chance to place limits on the size of such effects or to 



effectively parametrise them as for instance in Ref. [18|] . Moreover there are additional issues which may 



arise such as the role of small- j; effects that are neglected by conventional Qt resummations but may 
become important at values of x relevant for certain LHC processes such as Higgs production [l^ . Once 
again the existence or otherwise of such small- a; enhanced terms may be inferred from a comparison of 
conventional resummation with data in the region of large vector boson rapidity (which corresponds 
to small-x for one of the incoming partons). The need or otherwise, from a phenomenology viewpoint, 
of alternative theoretical constructions such as unintegrated parton densities [23]-[23| may also be 
discussed in the context of such a comparison. 

Further obscuring the issue of the relevant physics in the low Qt domain are certain technical 
issues concerning resummation of logarithms of M/Qt where M is the hard scale of the process. Most 
current methods rely on carrying out the resummation in impact parameter (6) space and then Fourier 
transforming to momentum space to obtain the Qt spectrum. Unfortunately the inversion to Qt space 
is not sensible unless one makes certain ad-hoc modifications to the pure resummed predictions. While 
these modifications are formally subleading from the point of view of the logarithms, they can become 
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numerically important in the Qt distribution. Clearly such artefacts do not pertain to genuine non- 
perturbative effects and hence, if significant, potentially obscure a meaningful comparison of theory 
and experiment. 

In summary the ever increasing accuracy of resummed predictions and the corresponding increased 
precision of the latest experimental data from hadron colliders as well as the need to use these techniques 
in the future (for instance in studying the distributions for new particles that may be discovered at the 
LHC) all contribute to the continued importance of studies of Qt spectra. 

In our present Letter we address the theoretical prediction, not for the Qt spectrum itself but 



consider instead the novel variables introduced in Refs. 23j, |2j] and recently measured by the DO 



collaboration 25|]. These variables labelled the ay and (j)* both crucially depend on the azimuthal angle 



between the final state leptons, at low Qt- Since it is possible to measure the qt and (j)* with 



significantly lower experimental errors 2j], in the low Qt region, an accurate theoretical prediction 



for these variables, along the line of those for Qt resummation, is important. The goal of our present 



Letter is to extend our existing resummed calculation for the qt distribution 26|] in such a manner as to 
provide a final matched result, taking into account resummation, including NNLL terms and matching 
it to fixed order predictions, accounting for the full DO experimental cuts. 

We organise the present Letter as follows. In the following section we introduce the (j)* variable 
measured by the DO collaboration and relate it in the soft limit to the qt variable treated in our earlier 



work 26|]. Having done so we provide in the subsequent section the resummed result for both the (p* 
and the ay, incorporating the cuts on the final state leptons as used by the DO collaboration. We also 
now include the a^L (NNLL) term in the resummed exponent and by expanding our results to order 



and comparing with fixed order predictions from the program MCFM 27|, demonstrate that we 
have complete analytical control of all large logarithms that may arise up to the two loop level. This 
fact allows us to straightforwardly match the fixed-order results from MCFM to the resummed results 
without the need to provide an ad-hoc prescription for treating subleading logarithmic pieces left over 
from the fixed order calculations. We conclude by pointing out that our matched resummed result can 
then be used as a basis for detailed phenomenology and comparisons between different theoretical and 
Monte Carlo approaches, which is work in progress. 

2. The 0* observable and its relationship to 

The (j)* observable as defined in Ref. [2J] is 

=tan((/>acop/2)sinr, (2.1) 

where i;^acop = tt — Ai;^ and Ai;^ is the azimuthal angle between the two leptons produced by the Z decay 
and sin 9* is the scattering angle of the dileptons with respect to the beam in the dilepton rest frame. 
Two variants of the (j)* arise by considering two different definitions of 9* where one can choose it to 
either be the angle in the Collins-Soper frame [2j (which is a particular dilepton rest frame) or one can 
simply boost along the beam direction such that the leptons make an angle 9* and 7^ — 9* with respect 
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to the beam |24]. This definition of 9* which avoids the necessity for measuring magnitudes of lepton 



momenta significantly helps the experimental resolution as explained in Ref. [24i |. From our viewpoint 
as far as the resummation is concerned, since we are interested in the low Qt region, where Qt is the 
pair transverse momentum, the various different definitions of 9* coincide. Let us consider the situation 
in the dilepton rest frame. In the limit Qt — > the angle 9* can be straightforwardly expressed (in 
terms of lab frame variables or boost invariant quantities) as sin^ 0* = 4/|,/M^, where It refers to the 
lepton transverse momentum and M is the dilepton invariant mass. 

Next we derive the relationship between the ut variable theoretically studied in Ref. 26|] and the 
0* defined above. We consider the variable 0acop/2 and write 




2 2 y 1 - cos A(/) ^ ^ 

Next one notes that 

Q2, = /2 ^ _^ /2 2 _^ 21tiIt2 cos A(j), (2.3) 

where Iti and It2 are the magnitude of the lepton transverse momenta. In the soft limit, Qt — >• 0, we 
have that Iti ~ It2 and in fact it is easy to see that in this limit the leptons are nearly back-to-back 
in the transverse plane which results in \Iti — It2\ = Qt \ cosa| where a is the angle made by the Qt 
vector with the lepton axis in the transverse plane. Using this in Eq. (12. 3p we get 



sin^ a « 2^ (1 + cos A0) . (2.4) 

where we have now set Iti ~ ^T2 = It- Also in the soft limit one can approximate 1 — cosAcf) ~ 2 so 
that we can eventually write 



^^acop_ / 1 + cos A(/> ^ Qt sin a 
2 ~ y 1 - cos A0 ^ 2/t ■ 

Recalling the definition of the ut variable as the component of Qt normal to an axis in the transverse 
plane which coincides with the lepton axis (in the Qt limit) one makes the identification qt = 
Qt sin a and hence we arrive at 

(/.*=tan(</-,,op/2)sinr (2.6) 

with M the dilepton invariant mass. 

Thus the (p* variable is straightforwardly related, in the small Qt limit to the ct variable resummed 
in our previous work. We can thus use our previous resummation and extend it to include the full 
experimental cuts for the (p* as well as the ut- Having done so we shall then check our resummation 
against fixed order predictions and confirm full analytical control of large logarithms up to the two- 
loop level, including the NNLL In 0* . Finally we provide the matching between the fixed-order 
corrections calculated up to two-loop level, and our resummation which puts our work into a form 
where phenomenological investigation of the non-perturbative piece can be directly carried out. 
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3. Resummed result 



Here we write down, without deriving it in detail, the resummed result which is a simple extension of 
the result we derived in our previous paper on the ut [2^ • The resummation is based on exponentiation 
of multiple soft and hard collinear radiation along with the factorisation of the observable's phase space. 
The latter point can be better expressed by noting that in the soft limit Qt — )• the (/>* (we avoid 
referring to the qt explicitly from now on) can be expressed in terms of the transverse momenta of 
multiple soft emissions as cp* ^ \ ^ ^ ^ sin cpi \ where (pi is the angle made by the soft emission ki with the 
lepton axis. Thus the condition for several soft gluon emissions to produce values of tan ((/>acop/2) sin 9* 
up to some fixed value 4>* , is simply (writing ku sin (pi = kyi ) 

0(^^*_^M^ =l^°°*,in(6Mr)ne^'"=-, (3.1) 

which shows the factorisation of the phase space constraint for the (p* in terms of the product of 
contributions from individual gluons. The above constraint on only a single component of the kT, 
kyi, is thus reflected by the presence of the sm(bAI(p*) function in contrast to the Bessel function that 
appears in the Qt resummation case, which shall account for the absence of a Sudakov peak in the (p* 
distribution, as explained in more detail later. 

a{(p*,M'^,cose*,y) =^ [ ^ sm{bM(P*) e'^^^'^^ {xi, X2,cose* ,bM, Hf) , (3.2) 

where 

M , - be'^'^ 

X12 = —^e^ and b = . (3-3) 

Vs 2 

The resummation is encapsulated by the b integral of the exponential suppression e"^^''^ with the 
sm{bM(p*) factor having arisen from the one-dimensional momentum conservation constraint as ex- 
plained above. We have written the result as differential in the dilepton invariant mass M as well as 
the dilepton rapidity (or equivalently the Z boson rapidity) y. In the following the factorisation scale 
Hf will be set equal to M. 

The function S has an identical structure to the Born level result, where 9* is the scattering angle 
in the dilepton rest frame. Note that S acquires a dependence on the impact parameter b because 
of the resummation of logarithms of b via DGLAP evolution, which then determines the scale of the 
parton distribution functions embedded in S. Including the contributions from the Z as well as from 
the virtual photon, we have 

S = (i + cos^9*){Ql-2QgViV,xi + {Af + Vi''){Al + v;^)X2)J'^ 

+ cos9*i-iQgAiAgXi + 8AiViAgVgX2)J'g , (3.4) 

where a sum over quark flavours q is implied. The above equation is naturally written as the sum 
of the two terms with different angular dependence: the first one is proportional to (1 + cos^0*) and 



4 



represents the parity conserving piece of the electro-weak interaction, while the term involving cos 9* 
is the parity violating piece. We notice that upon integration over the full 9* range, as well as over 
symmetric intervals, the parity violating term vanishes. The coefficients Ai^q and Vi^q are the electroweak 
couplings for lepton I and parton q, explicitly given by: 

Af=T^ and Vj = -2Qfs\n^ 9w. f = l,q, (3.5) 

where Tj? is the third component of the isospin. We also have introduced 

M2 (M2 - M|) 



Xi 



X2 = 



(M2 -M|)2 + r|M|' 



(M2-M|)2 + r|M|' 



(3.6) 



Ana 

In Eq. (|3.4p . J^^ are explicitly given by 

= {C0 h)q {xi,b) (C f2),-(x2,6) ± (C ® fi),-(xi,6) (C fa), (x2,6). (3.7) 

The convolutions involving the matrix of coefficient functions C and the vector of parton densities 
fi,2 for incoming hadrons 1 and 2 respectively can be explicitly written as 

(C fi), (x„ b) = 1^ ^Cq. (a. (0 , ^) ft (^z, i) , (3.8) 

where i = 1,2 and a sum over all flavours a is implied. 

The coefficient functions C represent perturbative corrections to the collinear branching of an in- 
coming parton a to a parton q which annihilates with q to form the Z boson. We note that the collinear 
enhanced terms generated by such a branching are incorporated to our accuracy into the scale of the 
pdfs fi via their dependence on the impact parameter b. Thus the coefficient functions represent only 
the non-logarithmic constant terms. For our purposes we shall only need the coefficient functions up 
to first order in a^: 

Cqaix, b) = 6q^6 (1 - x) + ^^^C^'J (x) + O (a^) . (3.9) 



The explicit expression for the coefficient functions is given in Appendix Appendix A 

We note that the argument of ag to be used while evaluating Ci should be of the order of the 
hard scale M if one works to NLL accuracy but is 1/6 if one wishes to achieve the NNLL accuracy we 
seek, hence we have retained an explicit b dependence of the coefficient functions. Other than the b 
dependence indicated in the coefficient functions and the pdfs the large logarithms in b are resummed 
into the exponential function involving R{b). 
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The resummation of logarithms in b space turns out to be precisely as for the case of Qt resumma- 
tion 12611 and hence the function R{b) is the same as for that case. R(b) can be expressed in the general 
form |28l ] 

R(b) = Lgi{asL) + g2{asL) + asg^{asL) ^ (3.10) 

with L = In (6^M^) . The functions gi and 52 describe leading and next-to-leading logarithms respec- 



tively while 53 refers to NNLL terms. A calculation by the Florence group [12l . Il5l ] extends up to the 
g^ term but a recent study by Becher and Neubert indicates that one of the assumptions involved in 
computing part of the (73 piece used in Qt resummation, namely that the three loop coefficient A'^^^ is 
identical to the one appearing in threshold resummation, may in fact not be valid [17,]- Since the piece 
of g'3 we refer to above starts at order a^, for the present Letter we shall avoid this issue pending further 
clarification by retaining only the well established a^L piece of g'3 which has been known for several 



years 29|]. We anticipate that the terms we neglect here corresponding to higher order contributions 
to g^ will in any case have a negligible effect, but note that it is trivial to include them if requirec|3- 



Explicit expressions for the functions appearing in the radiator are collected in Appendix A 

The DO collaboration recently measured the (j)* distribution with a set of cuts on the lepton momenta. 
More specifically, in the case of muons, they require the rapidity of both leptons to be |r7i^2| < 2 and 
transverse momenta pT of both leptons to be above 15 GeV. These cuts which leave the basic form 
of the resummation unchanged can be easily incorporated into our calculations by working out the 
appropriate limits for the 6* and y integrals and we have included them in our final result. 

In the following section we shall check the resummation by considering its expansion to order a1 and 



comparing with exact fixed-order results from MCFM 271] . We shall demonstrate that the resummation 
allows us analytical control over all logarithmic terms at a1 accuracy. 

4. Comparison to fixed-order calculations 

To expand the resummation to order a1 it is most convenient to consider the Mellin moments with 
respect to xi and X2. Then one can write the pdfs at scale 1/6 in terms of the pdfs at scale M using 
DGLAP evolution, which amounts to replacing the convolutions in Eqs. (j3.7p and (j3.8p by A^-space 
products as follows 

(C^fi) rxi,6)^|cfa, fiVAri%xp[-r(a„Af,)^]f(^i,M)l , (4.1) 



where C are the Mellin moments of the coefficient functions C, while Us in the argument of the two- loop 
anomalous dimension F is evaluated at scale M. We also need to expand (I/6): 



= a,(M) + pQal{M)L + O (af) . (4.2) 



similar conclusion was reached in a revised version of Ref. [l^, where it was found that the numerical impact of 
the modified jA*-^' is of the order of a few percent at small Qt- 
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The radiator R{b) can also be expanded up to the relevant in the form 



Rib) 



a 



2^ ,Gi2L2 + GuL) + (g;) (^23^=^ + G22L' + G21L) + O 



(4.3) 



where the coefficients Gnm are reported in Appendix A With the above information in place we can 
expand the result out to order straightforwardly and carry out the b integrals term by term to express 
the answer in </)* space. We also invert (in this case possible by inspection) the Mellin transform from 
Ni space back to Xi space. 

While the expressions for the expansion of the resummation up to order is too cumbersome to 
report here, we shall adopt the strategy of our previous work where we write down an expression for the 
difference between the integrated cp* and the very well known Qj- distributions. Because the logarithms 
in b space are the same in the two cases and because the sine and the Bessel function transforms have 
little effect on the expansion of the resummation, the logarithms in cp* mostly cancel the ones in Qt 
and the difference of the distributions admits a relatively simple analytical form: 



aiNuN2,(l>*] 



aiNi,N2,QT/2) 



Cro (A^1,A^2 



vr^d In^ 1 + ('-24C|C(3) - Svr^Cl - ^tt^Cf^o 



ro(iVi)fi(iVi)J f2,-(iV2) + 1 ^ 2 

fl„(iVl)f2„-(iV2) + 102 



In- 



(4.4) 



where only the first term involving In^ (^) had been checked against MCFM in Ref. [26] for the case of 
the aT distribution. In the above equation a{Ni,N2) represents the Mellin transformed cross-section, 
taking Mellin moments with respect to xi and X2 while ctq is the Born cross-section. Fq denotes the 
leading order anomalous dimension matrix and we have used the matrix notation we introduced earlier 
in Eq. (ISTH . 

By taking the derivative of Eq. ()4.4p . we can easily compute the difference between the logarithmic 
terms in the differential distributions: 



AD{e) 



1 



do d In e 



(^(0*) 



a {Qt/2) 



r/2=e 



(4.5) 



Subtracting AD from the corresponding fixed-order differential distribution D (0*) — D {Qt/2), com- 
puted with MCFM at NLO, we should find a result that tends to zero at small e. This difference is 
plotted in Fig. [1] and one notes that the result tends to zero as expected. Although we do not report 
them here, we have explicitly checked the expansion of the resummation against MCFM for the (/>*, the 
aT and the Qt distributions separately. 
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Figure 1: The difference between the NLO differential distributions for 0* and Qt/2 from MCFM after removal of 
logarithmic terms from the resummation. 



5. Matching to fixed-order 

Having demonstrated complete control over the divergent pieces at order we are in a position to 
propose a particularly simple matching formula for the differential (p* distribution, where one simply 
adds the resummation to MCFM results and subtracts the expansion of the resummation to order 
a^. This scheme is known as i2- matching [28i] and one does not need here to provide an arbitrary 
prescription to treat left over divergent pieces from the difference between MCFM and the expanded 
resummation. 

Before we can match the resummed result to fixed-order we have however to deal with the evaluation 
of the b integral in Eq. (13. 2p . As is well known, the radiator R(b), which resums terms logarithmic 
in b, is divergent both at small b and large b. The large b divergence is due to the Landau pole in 
the running coupling and is associated to non-perturbative behaviour. We shall hence impose a cut- 
off at a value 6max) chosen in the vicinity of the Landau pole, to avoid this problem. Additionally 
we need to handle the divergence at small b. Since the small b region describes large Qt or (j)* it 
is beyond the jurisdiction of our resummation formula. We decide to freeze the radiator R{b) for 
values of bM < 1, which is the region of large (f)* not controlled by the resummation [sO]- Different 
prescriptions, such as the one adopted by the Florence group, where one modifies the logarithms of 
b in R{b), In (JrM'^^ — )• In (l + IrM'^^ will be investigated in a future phenomenological study. Our 
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matched result corresponds to 



'i'?^* ) matched ^ d*/** ) resummed (d<?^*)MCFM ^'Z'* ) expanded 

where we have indicated by ( ^fe ) the expansion of the resummation to NLO which is achieved 

\ ^ J expanded 

by expanding the h space radiator to that accuracy and performing the h integral precisely as for the 
resummation. For the matching to be considered successful we ought to observe that the matched curve 
follows the resummed result at small (j)* while at large (j)* one expects the resummation to largely cancel 
against its expansion (up to relatively small terms varying as af) ^'^'^ hence the result should follow 
the NLO MCFM curve. 

Our final result for the matched differential (p* distribution is plotted in Fig. [21 together with the 
pure fixed order O (a^) calculation obtained from MCFM. We notice that the NLO calculation diverges 
in the region of small (f)* ^ while the resummed and matched results tend to a constant. This behaviour 
can be explained by performing a simple calculation at the double logarithmic accuracy: 



1 dcr _2£^1„2rf,2fl,f2^ 

o- d(j)* JO JO 



d(6M)cos(6Mr)e-"^'"'(^'*^')= / d{bM)e-^'^'(''^n 



1 + 



/*2 



vr 



F+0(<l)*^). (5.2) 



The above calculation indicates that the differential distribution goes to a constant at small (p*. In 
order to more accurately obtain the value of the asymptotic constant of the (p* distribution, one needs 
to include the full calculation up to at least NLL terms in the exponent, as we have done here, as well 
as include non-perturbative effects. 

While low values of (j)* or Qt can be obtained via Sudakov suppression or kinematical cancellation, 
in the present case of (p* the kinematical cancellation starts to dominate prior to the formation of the 
Sudakov peak, in contrast to the Qt case. This heralds the breakdown of the logarithmic hierarchy in 
the resummed exponent in the region corresponding to large values of b. Hence, a precise evaluation 
of the b integral becomes important, which however depends on a number of prescriptions referred to 
earlier. It would be interesting to study the implication of this observation for the phenomenology 
of non-perturbative effects and to compare the conclusions to the Qt case, as we intend to do in 
forthcoming work. 

From Fig. [2] we note that the matched curve tends to MCFM at large (p* but deviates from fixed- 
order significantly at lower (p* and follows the resummed curve as anticipated. We also note that the 
deviation from NLO will be significant over a large range in (p* where there are accurate data and 
hence we would see a visible role for both resummation and potentially non-perturbative effects in 
phenomenological studies. 
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Figure 2: The differential (j)* distribution computed at NLO with MCFM (dashed blue line) and our final result (solid 
red) obtained by matching our resummation, which contains some of the next-to next-to leading logarithms and hence is 
denoted by NNLL*. to the fixed order calculation at 0{o?s), from MCFM. The distributions are normalised to NLO cross 
section. 
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6. Conclusions and Outlook 



In this Letter we have presented a resummed result for the newly introduced 0* variable, closely 
related to the previously studied ax variable, which is accurately measured down to very low values of cp* 
by the DO collaboration. Our resummation is valid to full NLL accuracy but we also include the leading 
order a^L term of the NNLL resummation which simplifies the matching procedure to a large extent. 
We note that it is straightforward to include the full NNLL resummation and non-perturbative effects 
via a Gaussian smearing in b space which developments we leave to our forthcoming phenomenological 
study. 

For phenomenological purposes a number of further studies of the resummed distribution still require 
to be carried out. Primarily one needs to ascertain the stability of the resummation against the change 
of various parameters and prescriptions involved in generating the resummed and matched result. 
Amongst these parameters are just the traditional factorisation and renormalisation scales which here 
have been taken to be all equal and set to the invariant mass of the lepton pair. Additionally here one 
can explore the role of changing the resummation scale or in other words rescaling the resummation 
variable bM to assess the impact of subleading logarithms omitted by the resummation. However, with 
inclusion of full NNLL effects we may anticipate subleading logarithms to be essentially invisible. 

More importantly there are certain ad-hoc prescriptions inherent in b space resummation such as 
the treatment of the radiator at large and small b to avoid respectively the Landau pole and the 
spurious small b divergences beyond the control of resummation. For a complete phenomenological 
study therefore we intend to produce an uncertainty band on the resummed result which reflects the 
variation with respect to the change in these parameters and prescriptions. Once this is done we can 
focus on pinning down non-perturbative effects most accurately, for we which we shall also carefully 
explore the different rapidity (y) bins over which data exist. 

We believe that our accurate theoretical calculations will encourage the use of these novel variables, 
namely (p* and ut, for experimental studies of the Qt spectrum of the Z boson at the LHC. From 
a theoretical viewpoint it is worth noting that the study here is closely related to variables such as 



Acf) between jets measured recently by ATLAS [3l[ and CMS [33] and computed in Ref. [33|; those 
calculations therefore signify an extension of Qt resummation to processes with colour (jets) in the 
final state. Also for the future it will be interesting to compare and contrast different theoretical 
approaches including both the newly developed SCET approach [16] as well as Monte Carlo parton 
shower based predictions in the context of the cj)* variable and the DO data. 
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Appendix A. Explicit resummation formulae 

The radiator used in the resummation admits the following expansion: 

R{b) = Lgi{asL) + g2{asL) + aggsiagL), 



51 (A) 



Cf 



7r/?oA 



[-A-ln(l-A)] 



3*^^ ^ n ^ + (1 - A) - A)] CfPi 



53(A) 



27r^o 

Cl An^ - 48C(3) - 3 A 
TT^ 16 1- A 



27r2/32(i-A) 



A + ln(l-A) , 1^^2 



+ -ln^(l-A) 



+ 



(A.1) 

(A.2) 
, (A.3) 
(A.4) 



with A = as{M'^) PqL, L = Int^M^ and K = Ca - - frz/. We recall the reader that in this 

Letter we only keep the term in the NNLL coefficient gs which gives rise to single logarithms at a^. At 
the aimed accuracy we need to consider the running of the strong coupling constant at two-loops: 



l-as{M' 



,2,/3i In(l-p) 



^0 1-p J 



p = a,(M2)/3oln-^, 



where the coefficients of the QCD ^-function are defined as 



/3o 



IIC^ - 2n/ 



Pi 



The function R admits the following perturbative expansion: 

R{b) = g {GuL^ + GnL) + (g)' (Gsgi^' + G22L^ + G21L) , 
with the coefficients given by 



G12 
Gn 



3Cf 



G23 =C'F-7r^o 

G22 =CFiK-STrPo) 

G2i=C^(^7r2-12C(3)- 

+ CFnA---.^ 



+ CfCa ( 6C(3) - ^ + y 



(A.5) 

(A.6) 

(A.7) 

(A.8) 
(A.9) 

(A.IO) 
(A.ll) 

(A.12) 
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The explicit form of the coefficient functions is 

C,5(a,(l/6),x) =6{l -x) + ^^^Cf ((y " 4) '^(l - x) + {1 - x)^ + O (a^) , (A.13) 

C7,,(a,(l/6),x) =^i^x(l-x) + 0(af) . (A.14) 

(A.15) 

The matching to the fixed-order result computed with MCFM is most easily performed at the level 
of the differential distribution. Hence, the 6-integral we wish to perform are of the following form: 

Qj^-J -^s\n{bM(t)*)\u^ (buf = (t>*- j d{bM) cosibM cP*) In"" {bM)\ (A.16) 

Explicit results for m = 1,2,3,4 and (p* > are given by 



(i(6M)(/.*cos(6M(/.*)ln(6M)^ = -2, 



9 /"oo , 

-J d{bM)(t)* cos{bM(t)*)\a^ {buf = 41n(^4(/)* 

j d{bM)(l)* cos{bM(j)*)ln^ {buf = -Gln^ (^4(/)*2j - 27r^ 

^ d{bM)(t)* cos{bM(t)*) In^ {buf = Sln^ {^M*"^) + 8^^ In [acP*'^'^ + 128C(3). 



(A.17) 
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